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Abstract 

We calculate the isovector (flavor-nonsinglet) unpolarized quark- and antiquark distri- 
butions in the nucleon at a low normalization point in the large--/V c limit. The nucleon 
is described as a soliton of the effective chiral theory. The isovector distribution ap- 
pears in the next-to-leading order of the l/iV c -expansion. Numerical results for the 
quark- and antiquark distributions compare well with the parametrizations of the data 
at a low normalization point. This large-iV c approach gives a flavor asymmetry of the 
antiquark distribution (violation of the Gottfried sum rule) in good agreement with 
the measurements. 
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1 Introduction 



The parton distribution functions of the nucleon contain the non-perturbative information 
which enters in the cross section for deep-inelastic scattering and a variety of other hard 
processes. Their scale dependence in the asymptotic region is governed by perturbative 
QCD and well understood. The starting points of the perturbative evolution, however, that 
is, the distributions at a relatively low normalization point, belong to the domain of non- 
perturbative physics and can at present only be estimated using approximate methods to 
deal with the problem of the structure of the nucleon. 

A very useful approximation is the theoretical limit of a large number of colors, N c — > oo. 
It is known that in this limit QCD becomes equivalent to an effective theory of mesons, with 
baryons emerging as solitonic excitations At low energies the structure of the effective 
theory is determined by the spontaneous breaking of chiral symmetry. The first realization of 
the idea of the nucleon as a soliton of the pion field was proposed by Skyrme |2|, |3j , which was, 
however, based on an arbitrary choice of higher-derivative terms of chiral Lagrangian. A 
more realistic effective action for the pion field is given by the integral over quark fields with 
a dynamically generated mass, interacting with the pion field in a minimal chirally invariant 
way [|J. Such an effective action has been derived from the instanton vacuum of QCD, 
which provides a microscopic mechanism for the dynamical breaking of chiral symmetry pi. 
The chiral quark-soliton model of the nucleon based on this effective action || [7j has been 
very successful in describing hadronic observables such as the nucleon mass, iVA-splitting, 
electromagnetic and weak form factors etc. |§. 

Recently it has been demonstrated that it is possible to compute also the leading-twist 
parton distributions at a low normalization point in the chiral quark-soliton model of the nu- 



cleon H, [K], |TT] . This field-theoretic description of the nucleon allows to preserve all general 
requirements on parton distributions, such as positivity and the partonic sum rules which 
hold in QCD. In particular, it allows a consistent calculation of the antiquark distributions. 
The approach has by now been extended to transverse polarized distributions [12] as well as 



to off-forward quark distributions [13 



The large-iV c limit implies a classification of the parton distribution functions in "large" 
and "small" ones ||. Generally, at large N c the quark distributions are concentrated at 
values of x ~ 1/N C , and with the chiral quark-soliton model we aim to compute them for 
values of x of this parametric order. The distribution functions which appear in the leading 
order of the l/iV c -expansion are of the form 

D l ^ e (x) ~ N 2 C p(N c x), (1.1) 

where p(y) is a stable function in the large A^ c -limit, which depends on the particular dis- 
tribution considered. They are the isosinglet unpolarized and isovector (longitudinally or 
transverse polarized) distributions, which have been computed in Refs.[|, |K| [12|]. The isovec- 
tor unpolarized and isosinglet polarized distributions, on the other hand, appear only in the 
next-to-leading order of the l/iV c -expansion and are of the form 

D smal \x) ~ N cP (N c x). (1.2) 
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In this paper we compute the isovector unpolarized quark- and antiquark distributions 
in the chiral quark-soliton model. The general method for calculating the l/iV c -sub leading 
distributions and a discussion of their properties have been given in Ref . . Here we actually 
compute the distribution function, including the Gottfried sum, which is a measure of the 



flavor asymmetry of the antiquark distribution at the low scale [14, 15). This task requires, 



among other things, to generalize the analytical and numerical methods for the computation 



of distribution functions to the case of the l/iV c -sub leading distributions, taking into 
account the rotation of the chiral soliton. 

The isovector distribution considered here is of particular interest for understanding 
the relation of the low-scale parton distributions to the chiral effective dynamics. The 
distribution functions computed within the effective chiral theory correspond to distributions 
of "constituent" quarks and antiquarks, i.e., objects which possess a structure in terms of 
QCD quarks and gluonP] ||. The parameter governing the compositeness of the constituent 
quark is the ratio of the dynamical quark mass to the ultraviolet cutoff of the effective theory 
- the "size" of the constituent quark. In the instanton vacuum this ratio is proportional to 
the small packing fraction of instantons, (p/R) 2 . Before comparing the quark- and antiquark 
distributions computed in the effective chiral theory with the parametrizations of the data 
at a low normalization point |18| one should "resolve" the structure of the constituent quark. 
From the singlet distribution, u(x) + d(x) + u(x) + d(x), "resolving" the constituent quark 
structure one recovers the true singlet quark and gluon distributions. Phenomenologically 
one finds that gluons carry about 30% of the nucleon momentum at a normalization point 
of n ~ 600 MeV, so in the singlet case the "resolution" of the constituent quark structure 
is a rather sizable effect || |10J- In the isovector case considered here the change of the 



distribution due to the "resolution" of the constituent quark structure is expected to be less 
important, as this distribution does not mix with the gluon distribution and its normalization 
is scale independent (isospin sum rule). Consequently, the isovector distribution calculated 
in the effective chiral theory can almost directly be compared with the parametrizations of 
Ref. |TB[ , allowing one to draw conclusions about the model dynamics even in the absence of 



a complete understanding of the structure of the constituent quark. 

The plan of this paper is as follows. Section ^| contains a brief exposition of the effective 
chiral theory and the description of the nucleon as a chiral soliton, including the semiclassical 
quantization procedure. In Section |3] we derive the expressions for the isovector quark 
distribution functions in the effective chiral theory and discuss their properties. We first 
outline the expansion in the soliton angular velocity (l/iV c -expansion) which is necessary for 
computing the "small" (l/iV c -suppressed) distribution functions, in particular the isovector 
unpolarized one. We then discuss the important issue of ultraviolet regularization, which 



has been treated in detail in Ref.|T0[- We show that the distributions obtained in our 
approach satisfy the isospin sum rule. We also discuss the Gottfried sum. In Section | we 
briefly describe the numerical technique used for computation of the "small" 1 /^-suppressed 
distribution functions. Similarly to the "large" l/iV c -leading distributions they are computed 
as sums over quark single-particle levels in the background pion field; however, now one is 
dealing with double sums over levels. We then discuss the numerical results and compare 



1 For a more detailed discussion of this complex issue we refer to Refs. 16, 171; see also Ref. JIT 
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TBI. Conclusions and 



them to the parametrizations of the data at a low normalization point 
an outlook are given in Section [|. 

A calculation of the isovector unpolarized and isosinglet polarized distributions following 
the approach of Refs.[0, [TO] has recently been performed by Wakamatsu and Kubota |T9 



However, these authors have neglected certain contributions to the distribution function, 
which are important in particular at small values of x, as we shall discuss below. Also, a 
calculation of the isovector unpolarized structure function in a related approach has been 
reported in Ref. ||2U[|. In that calculation, however, only the contribution of the so-called 
valence level is taken into account. This approximation leads to a number of inconsistencies, 
as has been discussed in Refs. 0. ITC 



2 The nucleon as a chiral soliton 

The starting point for the chiral quark-soliton model of the nucleon is the effective action 
for the pion field, which is obtained by integrating over quark fields in the background pion 

field ma. 



exp(iS&[U(x)]) = J DiPD^j exp i J d 4 x^(i$ - iW77 75 )</> 



(2.1) 



Here, ip is the fermion field, M the dynamical quark mass, which is due to the spontaneous 
breaking of chiral symmetry, and the pion (Goldstone boson) field is described by an SU (2) 
matrix, U(x), with 

U^(x) = l4pf/(x) + ^pV( x ). (2.2) 

In the long- wavelength limit, expanding in derivatives of the pion field, the effective action 
Eq. (|2.1| ) reproduces the Gasser-Leutwyler Lagrangian with correct coefficients, including the 
Wess-Zumino term. It is understood that the effective theory defined by Eq.( |2.1| ) is valid for 
momenta up to an UV cutoff, which is the scale at which the dynamical quark mass drops 
to zero. We shall take in the discussion here the quark mass to be momentum-independent 
and assume divergent quantities to be made finite by applying some UV regularization later. 
[Why this is generally justified will be discussed below.] 

The effective action Eq. ( |2.1| ) has been derived from the instanton vacuum, which provides 
a natural mechanism of dynamical chiral symmetry breaking and enables one to express the 
parameters entering in Eq. ( |2.1| ) — the dynamical mass, M, and the ultraviolet cutoff — in 



terms of the QCD scale parameter, Aq C d @- In particular, the cutoff is given by the average 
instanton size, p _1 ~ 600 MeV 

In the effective chiral theory defined by Eq. ( |2.1|) the nucleon is in the large- N c limit 
characterized by a classical pion field ( "soliton" ) . In the nucleon rest frame it is of "hedgehog" 
form 0, 



Z7 c (x) = exp [z(n • r)P(r)} , 



where P(r) is called the profile function, with P(0) = —tt and P(r) — * for r — > oo. Quarks 
are described by one-particle wave functions, which are determined as the solutions of the 
Dirac equation in the background pion field, 



H(U c )\n) = E n \n). 



(2.4) 



Here H(U C ) is the single-particle Dirac Hamiltonian in the background pion field given by 
Eq.(P), 



H(U) 



(2.5) 



The spectrum of H(U C ) includes a discrete bound-state level, whose energy is denoted by 
E\ey, as well as the positive and negative Dirac continuum, polarized by the presence of the 
pion field. The soliton profile, P(r), is determined by minimizing the static energy of the 
pion field, which is given by the sum of the energy of the bound-state level and the aggregate 
energy of the negative Dirac continuum, the energy of the free Dirac continuum (U = 1) 
subtracted [HI, 



tot 



Nr. 



E E n - £ 



occup. 



occup. 



N c E lcv + N c ( E n -E^), 

n 

neg.cont. 



(2.6) 



and in the leading order of the 1 /^-expansion the nucleon mass is given simply by the value 
of the energy at the minimum, 



M 



N 



E 



tot 



(2.7) 



u=u c 



The expression for the energy of the pion field, Eq.( [2.6| ), contains a logarithmic ultraviolet 
divergence due to the contribution of the Dirac continuum and requires regularization. In the 
calculation of the parton distributions below we shall use a Pauli-Villars regularization (see 



Subsection |3.2|) , so we regularize also the energy using the Pauli-Villars method. Following 
Ref.iTT! we define 



E 



tot,reg 



N c E lev + N c 



E (En-E^) 



n 

neg.cont. 



M 2 



E (e 



PV,n 



E 



(0) \ 
PV,n) 



n 

neg.cont. 



, (2.8) 



where Epy^ n are the eigenvalues of the Hamiltonian, Eq.( |2.5| ), with M replaced by a regulator 
mass, Mpv, and the same pion field. The value of the regulator mass, which now plays the 
role of the physical cutoff, can be fixed from the pion decay constant, 



N C M 2 

471-2 



log 



Mlv 
M 2 ' 



(2.9) 



6 



[Numerically, F n = 93 MeV.] Note that we do not subtract the finite contribution of the 
discrete level for the spectrum with the Pauli-Villars mass; this prescription leads to a stable 
minimum of the regularized energy functional, Eq. fl2.8Q , with respect to the profile function 
[ II] . The soliton profile (the "self-consistent" pion field) for this UV regularization has been 



determined in Refs. [|TT|, (2l|[ . 

In higher order of the l/iV c -expansion one must take into account the quantum fluctua- 
tions about the saddle-point pion field. A special role play the zero modes of the pion field. 
In fact, the minimum of the energy, Eq. fl2.6p , is degenerate with respect to translations of 
the soliton field in space and to rotations in ordinary and isospin space. For the hedgehog 
field, Eq.( |2.3|) , the two rotations are equivalent. Quantizing the zero modes modes, i.e., in- 
tegrating over collective rotations and translations, gives rise to nucleon states with definite 
momentum and spin/isospin quantum numbers ||, ||. One performs a (time-dependent) 
rotation of the hedgehog field, Eq.(|2.3j), 



C/ C (x) -> R(t)U c (x)R\t), (2.10) 

where the collective coordinate, R(t), is a rotation matrix in S'L r (2)-flavor space. The func- 
tional integral over the collective coordinate can be computed systematically within the 
l/iV c -expansion. The moment of inertia of the soliton is 0(N C ) (the nucleon is "heavy"), 
hence the angular velocity is 

Q = fi a y = -iR}R = oQr) (2.11) 

(the collective motion is "slow"), and one can expand in powers of the angular velocity. To 
leading order in Q the collective motion is described by a Hamiltonian 

H TOt = — — = —2-, (2.12) 

rot 2/2/ V ^ 

where S a and T a are the right and left angular momenta, and the Hamiltonian Eq. ( |2.12| ) has 
been obtained by the "quantization rule" 

fi Q - y- (2.13) 

Here / denotes the moment of inertia of the soliton. It can be expressed as a double sum 
over quark single-particle levels in the background pion field, 

N c ^ ^ (n\T a \m) (m\T a \n) 

i = -i E £ - ' _ ' ' - ■ (2-14) 

occup. non— occup. 

Here the sum over n runs over all occupied states, i.e., the discrete level and the nega- 
tive Dirac continuum, the sum over m over all non-occupied states, i.e., the positive Dirac 
continuum. 

The Hamiltonian Eq. fl2.12 ) describes a spherical top in spin/isospin space, subject to the 



constraint S 2 = T 2 , which is a consequence of the "hedgehog" symmetry of the static pion 
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field, Eg. (|2.3D. Its eigenf unctions, classified by S 2 = T 2 , S3 and T 3 are given by the Wigner 
finite-rotation matrices ||, 

<P S S S(R) = V2STT(-ir+ T -W s ^ S:j (R). (2.15) 

The four nucleon states have S = T = 1/2, with S3, T3 = ±1/2, while for S = T = 3/2 
one obtains the 16 states of the A resonance. The rotational energy, S(S + 1)/(2J), gives a 
l/iV c -correction to the nucleon mass, which should be added to Eg. fl2.7l ). m particular, the 
nucleon-A mass splitting is given by 

M A -M N = A (2.16) 



The expression for the moment of inertia, Eq. (|2.14 ), contains an ultraviolet divergence 



which is to be removed by the ultraviolet cutoff. The ultraviolet regularization of the moment 
of inertia must be consistent with that of the isovector guark distribution function; it will 
be discussed below in Subsection W71. 



3 The isovector unpolarized quark distribution 

3.1 The isovector distribution function in the large— N c limit 

To compute the twist-2 guark and antiguark distribution functions in the effective chiral 
theory we start from their "field-theoretic" definition as forward matrix elements of certain 
light-ray operators in the nucleon, which can be regarded as generating functions for the 
local twist-2 operators |p2| . Alternatively, one could start from the "parton model" definition 
as numbers of particles carrying a given fraction of the nucleon momentum in the infinite- 
momentum frame P3| - both ways lead to identical expressions for the guark distribution 



functions in the chiral guark-soliton model ||, |10|- The unpolarized distribution (/ denotes 
the guark flavor) is given by 



D nnpolJ (x) = i- Jdz-e^-(P|^(0)7 + ^/W|P> 



2+=0, 2 ± =0 

£> nnpo i,f(x) = -{x^-x}. (3.1) 
Here, z ± and ^ denote the usual light-like vector components and Dirac matrices, 

± z" ± z* ± V±7 3 

2 = ~VT' 7 = —JT- (3 ' 2) 

In the longitudinally polarized case one should replace in Eq.( [3.1|) 7 + — > 7 + 75, and the 
polarized antiguark distribution is given by the function at x — > — x without minus sign. In 
Eq . (|3~T| ) we have dropped the gauge field degrees of freedom (the path-ordered exponential 
of the gauge field), a step which is justified when working in leading order in the ratio of M 
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to the UV cutoff, viz. the packing fraction of the instanton medium; see Refs. |Tfi|, |T7 for a 
detailed discussion. 

Thanks to its relativistically invariant definition the matrix element Eq. (|3.1| ) can be 
evaluated in any frame; for us it is convenient to compute it in the nucleon rest frame. The 
calculation follows the usual procedure for computing matrix elements of quark bilinears 
in the chiral quark-soliton model within the 1/N C expansion. The limit N c — > oo justifies 
the use of the saddle point approximation for the pion field. The matrix element can be 
calculated with the help of the quark Feynman Green function in the background pion field, 



G F (y°, y] x ,x) 



(y°,y\[id t -H(U))- 



\x 



,x>. 



Here the saddle-point pion field is the slowly rotating hedgehog field, Eq. (|2.10|) 
ansatz the Green function, Eq. ( |3.3j ), takes the form 

[idt ~ H(U)}- 1 = R(t) [id t ~ H{U C ) - n]- 1 R\t), 



(3.3) 



For this 



(3.4) 



where Q is the angular velocity, Eq. (|2.11 ). The matrix element of a quark bilinear between 
nucleon states of given spin and flavor quantum numbers is obtained by integrating over 
soliton rotations with wave functions of the collective coordinates, R, corresponding to a 
given spin/isospin state. In addition, one has to perform a shift of the center of the soliton 
and integrate over it with plane-wave wave functions in order to obtain a nucleon state of 
definite three-momentum. This leads to the following expression for the forward matrix 
element of a color-singlet quark bilinear in the nucleon state (T denotes the time-ordered 
product): 

(P = o, S = T, S 3 , T 3 \ T {^{x)T^(y)} |P = 0, S = T, S 3 , T 3 > 

R(T)=R 2 

= -2iM N N c [ dRx [ dR 2 fef (ifc " 



SfeORi) / DRDet[idt-H{U e )-n] 



x / d 3 X Tr 



R(-T)=R 1 

R\x°)TR(y°) (y°, y - X| [id t - H(U C ) - n}' 1 \x°, x - X } 



(3.5) 



Here T denotes a matrix in Dirac spinor and isospin space, and Tr . . . implies the trace over 
Dirac and flavor indices (the sum over color indices has already been performed). The path 
integral over R(t) can be computed using the fact that in the large iV c limit the angular 
velocity of the soliton, Q = —iR^R, Eq. (|2.11|) , is suppressed, which allows to expand both 



the Dirac determinant and the propagator in the integrand in local powers of derivatives of 
R(t). In particular this expansion gives rise to the kinetic term 



exp 



(3.6) 



where I is the moment of inertia of the soliton, Eq.( j2.14 ). The path integral over R(t) with 
this action can now be computed exactly; it corresponds to the rigid rotator described by 
the Hamiltonian Eq. ( |2.12[) . 
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When expanding the integrand of the path integral in Eq.( |3.5| ) in powers of f2 we keep all 
linear terms in Q in addition to the exponentiated kinetic term. The calculation of the path 
integral over R(t) is then equivalent to replacing the angular velocity by the spin operator, 
S, according to the "quantization rule", Eq. Q2.13l ). When dealing with the matrix element 
of a non-local bilinear operator as in Eq.( [3.5|) , terms linear in Q in the integrand of the path 
integral arise from two sources: 



Expansion of the quark propagator (cf. Eq. (|3.4j )): 

[id t - h(u c ) - n] 

= [id t - H(U C )}- 1 + [idt - H{U c )\~ l Q [id t - HiU,)}- 1 + ... 



(3.7) 



Expansion of the nonlocal object R j; (x )TR(y ). This expansion can be performed in 
two alternative ways: 



or 



R\x )TR{y°) 
tf(x°)TR{y°) 



rf(x°)TR(x°) + (y° - x°)R^x )TR(x°) + ... (3.8) 



R\y°)TR(y ) - (y° - x°)R\y°)TR(y ) + 



(3.9) 



It can be shown [0] that both choices lead to the same result for the rotational correction 
to the matrix element. 

Turning now to the calculation of parton distribution functions, defined as matrix ele- 
ments of the non-local operator Eq^O), it was shown in Ref. || that the isosinglet unpolar- 
ized distribution, u(x) +d(x), is non-zero already in the leading order of the l/iV c -expansion, 
i.e., at order Q° in the expansion of the integrand of Eq.( [3.5|) . The isovector unpolarized dis- 
tribution, on the other hand, appears only in the next-to-leading order, after expanding the 
integrand to order Q 1 . It consists of two pieces, which arise, respectively, from the expansion 
of the propagator, Eq. fl3.7p, and due to the non- locality of the operator, cf. Eq.(|3.8|). 

u(x)-d(x) = [u(x)-d(x)]M + [u(x)-d(x)f\ (3.10) 

The first contribution from the expansion of the propagator, Eq.( [3.7|) , is given by 

iM N N c 



[u(x) -d(x)} {1) 



s 3 . 



X (z 







dz"e 



XI 



„ixMpfZ 



Air 

'° J d 3 X J dR 



h T=S 



(R) 



Tr 



id t ~ H(U C 



_ S a a 
21 



id t - H(U C 



R^t 3 R (1 + 7 V) 



|0,-X) 



(3.11) 



The second contribution, originating from the non-locality of the operator, Eq.( |3.8| ), is given 
by 



[u{x)-d{x)f ] 



iM N N c 



Ait 
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oo 

$3 —oo 



X (z°,z-X\ 



{] < " M » z0 I d 3 x 1 (in 
1 



lT=S 



Tr 



21 



z°R^r 3 R —S a r a (1 + 



7V) 



id t - H(U C 



|0,-X) 



,z°=— z" z 1 



(3.12) 



Here S" a is the spin operator, acting on the rotational wave functions, which enters through 
the "quantization rule", Eq. fl2.13| ). Now it is a matter of routine calculations, typical for the 
chiral quark-soliton model, to compute the rotational matrix elements in Eqs.( [3.11 , p. 12 



and evaluate the functional traces using the basis of eigenfunctions of the Dirac Hamiltonian 
Eq. (|2.5| ). The details can be found in Appendix |A|. The result for the first contribution 
Eq.(|£IH), is 

N r M N 



[u(x) - d(x)] {1) = (2T 3 )- 



12/ 



x 



2 E E 



(n\T a \m)(m\T a (l + -f°-f 3 )S(E n + P 3 - xM N )\n) 



n m J - , m 

OCCUp. Em^E n 



E m — E r , 



E (n\T a \m)(m\T a (l + 1 ° 1 3 )5\E n + P 3 -xM N )\n) 



(3.13) 



occup. E m =E n 



where T 3 = ±1/2 is the isospin projection of the nucleon. Eq.( 



expresses the isovector 



distribution as a sum over quark single-particle levels in the background pion field, cf. 
Eq.(|J). In the first term on the R.H.S. of Eq. flST^) the outer sum over n runs over all 
occupied levels, including the bound-state level as well as the negative Dirac continuum, 
while the inner sum runs over all levels with the restriction that E m 7^ E n . 2 This formula 
assumes a quasi-discrete spectrum and can be directly applied to numerical calculations in 
a finite box; for the continuum case see Appendix The second contribution, Eq. (|3.12|) , 
can be expressed through the derivative of the isosinglet quark distribution, calculated in 
leading order of the l/iV c -expansion: 



[u(x) -d(x)] {2) 



-(2T 3 )^x^ J2 (n\(l + j°j 3 )6(E n + P 3 -xM N )\n) 



occup. 



•(2T 3 



w x ^ Wx)+rf(x)] 



leading 



(3.14) 



We note that, as in the case of the isosinglet distribution P, PH|, there exists an equivalent 
representation of the isovector distribution as a sum over non- occupied levels (see Appendix 



[u(x) -d(x)] {1) 



-(2T 3 



N r .M 



A' 



12/ 



2 Here, as in Eqs.(2.6, 2.14j ), the sums over n, m denote the sum over all quantum numbers characterizing 
the single-particle states; note that the energy eigenvalues are in general degenerate in the third component 
of the "grand spin" , the sum of the quark angular momentum and isospin H . 
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X 



2 E E 



(n\T a \m)(m\r a (l + -f°-f 3 )6(E n + P 3 - xM N )\n) 



E — E 

n m J -'m J -'n 

non— occup. E m j^E n 



"EE (n\T a \m)(m\T a (l + ^ 3 )5'(E n + P 3 - xM N )\n) 



non— occup. E m =E n 



(3.15) 



[u{x) -d(x)f ] 



(2T3)^ x #- 



E 

n 

non— occup. 



(n\(l + 7°7 3 )<K£n + P 3 ~ xM N )\n). 



(3.16) 



In both Eqs. (|3.13| , p. 14|) and ( |3.15| , |3.16| ) vacuum subtraction, i.e., subtraction of the 
corresponding sums over vacuum levels (U = 1) and with vacuum occupation numbers, is 
understood. It can be shown that vacuum subtraction is required only for x < in the sum 
over occupied, and for x > in the sum over non-occupied levels. This fact is important for 
the numerical calculations: we shall use the representation as a sum over occupied states for 
x > 0, and as a sum over non-occupied states for x < (see below). 

The two contributions to the isovector distribution function, Eq. (|3.13|) and Eq. (|3.14j) , 



which have emerged from the expansion of the integrand of Eq.([J3), both contain delta 
function type singularities at x — 0. For the second contribution, Eq. ( |3.14j ), this is immedi- 
ately obvious: the isosinglet distribution is discontinuous at x — 0, so its derivative contains 
a delta function centered at x — 0. In a similar way one can convince oneself that also 
the second term in Eq. (|3.13|) exhibits such a singularity. These delta function singularities 
cancel when the two contributions are added, cf. Eq. ( |3.10| ). For the numerical calculations 
it is convenient to regroup the terms in such a way that this cancellation happens at the 
level of analytical expressions, i.e., that the singular terms are combined in one term. By 
inserting intermediate states this term can be expressed as a double sum over levels, similar 
to the first term on the L.H.S. of Eq.( |3.13D . Doing this simple rearrangement one obtains 



[u{x)-d{x)\ = [u(x) - d(x)] (V) + [u(x) - d(x)] (2 '\ 



(3.17) 



[u{x)-d{x)f' ) = (2T 3 ) 



3n N c M n 



x £ E 

n m J -'m J -'n 



occup. E m ^E. 

[u{x] -d(x)} {2,) 
d 



-(2T A ) 



61 

(n\T a \m)(m\T a (l + 7° 7 3 )5(£n + P 3 - xM N )\n), (3.18) 



121 



x -pT E E (n\T a \m}(m\T a (l + 1 ° 1 3 )5(E n + P 3 -xM N )\n), 
ox n m 

OCCUp. Em^E n 



(3.19) 



and a corresponding representation as sums over non-occupied states, cf. Eqs.( p.l5| , |3.16| ). 
The expressions for the distribution function Eqs. (|3.17| , [3.18||3.19] ) will be used in the actual 
numerical calculations. 



12 



We would like to comment on a recent calculation of the isovector distribution in the 
chiral quark-soliton model by Wakamatsu and Kubota [[HJ, who follow basically the same 
approach as we do here. However, in the expansion of the integrand of the path integral 
over soliton rotations, Eq.(|3.5), these authors drop the terms arising due to the non-locality 
of the quark bilinear operator, Eq.( 3.19 ). They argue that the non-locality of this operator 
is of the order of the inverse large momentum transfer in deep-inelastic scattering, 1/Q, so 
that the slow rotational motion of the soliton does not affect the operator. This argument, 
however, seems not appropriate: the light-like separation between the quark fields in the 
operator for the parton distribution is governed by the Bjorken variable x (more precisely, 
by xMn), not by the momentum transfer Q 2 . After factorization of the DIS cross section in 
QCD the hard momentum enters only in the coefficient functions, which can be computed 
perturbatively, not in the operator for the parton distribution functions. Hence there is no 
reason to drop the contribution Eq.( |3.19| ) from the Q-expansion, the more that the non- 
locality of the quark bilinear operator is taken into account at other places in the calculation 
of the parton distribution. 

The authors of Ref.|19| also argue that in the contribution Eq. ( |3.18|) the double sum 
over levels should be restricted to include only transitions from occupied to non-occupied 
states, since transitions from occupied to occupied states, as are present in Eq. ( |3.18| ) (cf. the 
derivation in Appendix [A]) would violate the Pauli principle. While this argument is correct 
for local operators, we are dealing here with an operator non-local in time, for which this 
so-called Pauli blocking is not restrictive. We thus see no reason for dropping the occupied- 
to-occupied transitions in the double sum over levels in Eq. (|3.18|) . Indeed, it turns out that 
the full expression Eq.( [3.18j ) is required to ensure equivalence of summation over occupied 
and non-occupied states, as will be discussed below. 

The dropping of the contribution Eq.( |3.19D to the parton distribution function in Ref.|T9 



has rather drastic consequences. One may expect the differences to the full result to be most 
significant at small values of x, since they correspond to large light-like separations of the 
bilinear operator. To illustrate the differences we have performed a numerical calculation of 
the distribution function following the prescription of Ref . [[lSfl and compared the results to 
those of our calculation, where we take into account all terms. [The details of the numerical 
calculations and the ultraviolet regularization are given below in Subsection |3.2| and Section 
|].] Fig.|I] shows the contribution of the negative Dirac continuum to the sum over n in in 
Eqs.( |3.18| , |3.19| ) in both cases. One sees that the two distributions differ strongly. The full 
result (solid line), corresponding to the sum of Eqs.( |3.18"D and ( |3.19|) , is large at x = and 
vanishes rapidly for larger values of x. The result obtained after dropping Eq.( |3.19| ), as done 
in Ref. |19|, is shown by the dashed line; in this case the distribution is small at x = and 



does not vanish for large values of x. 

To conclude this discussion, we note that the truncations made in Ref. |19| seem problem- 
atic in yet another respect. The calculation of the parton distribution in Ref. [[19] starts from 
the bilinear operator ip . . .ip\ however, one may just as well take the opposite ordering of 
the operator, —ip ...-?/>. In QCD the two orderings give equivalent expressions for the parton 
distribution, thanks to the anticommutativity of the fermion fields at space-like separations. 
In the chiral quark-soliton model, the two orderings lead, respectively, to representations 
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of the parton distribution as sums over occupied and non-occupied quark single-particle 
states. The equivalence between these two representations holds only for the full expression 
for the distribution as it derives from Eq. (|3.5| ), i.e., for the sum of all terms in the expansion 
with respect to angular velocity; see Appendix Consequently, dropping the contribution 
Eq.( 3.19|) (or otherwise truncating the sums over levels) one violates this equivalence. We 
have verified this in the numerical calculations. For the full result (the solid line in Fig.[TJ) 
we observe equivalence of summation over positive and negative energy states, but not for 
the distribution computed according to the prescription of Ref. [Ll| (the dashed line in Fig.|TJ) 
when computed by summing instead over positive-energy states. 



3.2 Ultraviolet regular izat ion 

In the previous section we have obtained expressions for the isovector quark distribution 
as sums over quark single particle levels in the background pion field. In Eqs.( |3.18t |3.19| ) 
the sum over n runs over all occupied levels, that is, the discrete bound-state level and the 
negative Dirac continuum, schematically 



E 

n 

occup. 



lev 



+ E 



(3.20) 



n 

neg.cont. 



where the ellipsis denotes the "inner" sum over m with the restriction that E m ^ E n . We 
shall refer to the two terms on the R.H.S. of Eq.( |3.20| ) as the "level" and "continuum" 
contributions. One should keep in mind, however, that this distinction is quite formal, and 
that only the total sum over occupied (or, equivalently, non-occupied) states has physical 
significance. 

In the derivation of the expressions for the quark distribution function we have so far 
not taken into account the ultraviolet cutoff intrinsic in the effective chiral theory. In fact, 
the expressions for the isovector distribution, Eqs. (|3.13| , |3.14|) viz. Eqs. (|3.18| , |3.19|) , contain 
an ultraviolet divergence due to the Dirac continuum contribution, and thus require regu- 
larization. Quite generally, in a calculation of parton distribution functions there are very 
strong restrictions on how one should introduce the UV cutoff in the effective theory. The 
point is that one has to preserve certain general properties of parton distributions such as 
positivity, sum rules etc., which can easily be violated by an arbitrary UV regularization. 
Specifically, the regularization should preserve the completeness of the set of quark single- 
particle wave functions in the soliton pion field. One possible regularization method which 
fulfills all requirements is a Pauli-Villars subtraction, which was used in the calculations 
of the iV c -leading distributions in Refs.|9|, [10], p| . We shall also employ this method here. 
Thus, we regularize Eqs.( |3.13] , |3.14j ) viz. Eqs.( |3.l8| , |3.19j ) as follows: 



\u[x) 



d(x)\ 



reg. 



In = lev,Af 



E 

n 

neg.cont. 



M 2 
M 2 PV 



E 

n 

neg.cont. 



L" • -IM PV ' 



(3.21) 



where the subscript M PV denotes the corresponding expression computed with the con- 
stituent quark mass, M, replaced by the regulator mass, M PV , cf. Eqs.(|2~8|, |2T^ ). This 
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subtraction removes the logarithmic divergence of the distribution function, as can be shown 
using methods similar to those developed in Ref.[T(| (gradient expansion). We do not in- 
clude the contribution of the discrete level to the isovector distribution in the Pauli-Villars 
subtraction; this contribution is finite and does not need to be subtracted. 



The moment of inertia of the soliton, Eq.( 2.14 ), is also ultraviolet divergent. It must be 



regularized consistently with the isovector distribution in order to preserve the isospin sum 
rule (see Subsection |3l|). This is achieved by regularizing it by a Pauli-Villars subtraction 
analogous to Eq. ( |3.21| ). 

We stress again that the calculation of parton distributions in the effective chiral theory, 
i.e., the identification of the twist-2 QCD operator with an operator expressed in terms 
of fields of the effective theory, Eq.( |3.1|) , is based on the parametric smallness of the ratio 
of the constituent quark mass, M, to the ultraviolet cutoff (here: Mpy). The approach 
is consistent in the sense that the expressions for the parton distributions in the effective 
theory are logarithmically divergent and thus insensitive to the details of the ultraviolet 
regularization (assuming that the latter does not violate any important properties such as 
completeness etc.). It should be noted, however, that the lack of precise knowledge of the 
"true" ultraviolet regularization leads to a theoretical uncertainty at the level of finite terms, 
oc M 2 /Mp V . For instance, Pauli-Villars subtraction of the finite level contribution would 
correspond to a modification of the distribution function by an amount of order M 2 /Mp v . 
The numerical results presented in Section ^ should be understood as a calculation of the 
parton distributions with a typical regularization fulfilling all general requirements, with the 
theoretical uncertainty in the ultraviolet regularization leading to a numerical uncertainty 
of typically 10-20 % — except for the sum rules, which follow from general principles not 
violated by the ultraviolet regularization. 

3.3 The isospin sum rule 

We now want to demonstrate that the expressions derived for the isovector distribution in 
the chiral quark-soliton model are consistent with the isospin sum rule, 

1 i 

J dx[u(x) - d{x)\ = J dx[u{x) - u(x) - d(x) + d(x)] = 2T 3 . (3.22) 

-1 o 

Consider the representation of the isovector distribution function as a sum over quark levels, 
Eqs.( p.l8| , |3.19| ). Since the contribution [u(x) — d{x)}^ 2 '\ Eq. ( |3.19| ), is a total derivative, it 



does not contribute to the integral in Eq.( 3.22|) , and we can concentrate on the contribution 
[u(x) — d(x)}^ 1 >, Eq.( |3.18| ). Integrating Eq.( |3.18| ) over x, replacing in the large-iV c limit the 
integral from —1 to 1 by the integral over the whole real axis, we use up the delta functions 
in Eq. fl3.19p and obtain 



, r / \ ,/ m / ^i\N c ^ ^ (n\T a \m) (m\T a \n) 
dx[u(x) - d(x)] = (2T 3 )-^ J2 E I _ J, ■ (3-23) 

_1 OCCUp. Em^En 
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The hedgehog symmetry of the pion field, Eq.( |2.3| ), has allowed us to drop the 7°7 3 -term. 
Taking into account that the soliton moment of inertia, Eq. ( 12.14]) , can equivalently be written 

as 

I = ^£ J2 ^ ra l m ^ m l ra ^ (3 24) 

OCCUp. -C'm 

we immediately reproduce the isospin sum rule, Eq.( |3.22| ). Also, it can be seen that the 
Pauli-Villars regularization of the distribution function and the moment of inertia as defined 
by Eq. (|3.21|) does not upset this proof; in particular, it preserves the equivalence of the two 
representations of the moment of inertia, Eq.( [2.14| ) and Eq. ( p. 24] ) . 



3.4 The Gottfried sum 

An interesting quantity related to the isovector antiquark distribution is the Gottfried sum 



14], [L5|] , which is defined as 

Ig = I + I f 1 dx [u(x) - d(x)\. (3.25) 
3 3 Jo 

The integral on the R.H.S. is scale-dependent only at two-loop level; its scale dependence 
is negligible over the entire perturbative region. If the sea quark distribution were isospin 
symmetric, u(x) = d(x), which would be the case if, for example, one assumed the sea 
quark distribution to be generated entirely radiatively, this quantity would be equal to 1/3 
(Gottfried sum rule). However, the NMC experiment [^4] finds a significant deviation from 
this value, 

I G = 0.235 ±0.026 at Q 2 = 4GeV 2 , (3.26) 

indicating that the sea quark distribution is rather far from flavor-symmetric. Note that 
the Gottfried sum rule does not follow from any fundamental principles of QCD. In fact, 
the large- N c picture of the nucleon as a chiral soliton naturally explains the presence of a 
flavor-nonsymmetric antiquark distribution. 

The expression for the Gottfried sum in the chiral quark-soliton model is obtained by 
integrating the expression for the total isovector quark distribution, Eqs.( |3.18| , |3.19| ), over 
x from —1 to 0, keeping in mind that the antiquark distribution is given by minus the 
expressions in Eqs.( |3.18| , |3.19| ) at — x. Since the integration extends only over half of the x- 
axis, also the total derivative term, Eq.( |3.19|) , contributes to the Gottfried sum, in contrast 



to the isospin sum rule, Eq.( |3.22| ), where the integration runs over the whole x-axis and 
Eq. (|3.14|) drops out. This contribution to the Gottfried sum in the chiral quark-soliton 
model has not been mentioned in Ref. 0. 

We refrain from writing down the analytic expressions for the Gottfried sum obtained by 
integrating Eqs.( p.l3"j |3.14|) . Instead, we shall compute Eq.( [3.25|) by integration of the nu- 



merically computed antiquark distribution function, see below. In particular, the ultraviolet 
regularization of the Gottfried sum and all other properties will follow directly from that of 
the distribution functions. 
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We note that the Gottfried sum in the large-iV c limit has been discussed previously in 

These studies 



the context of the chiral quark-soliton model |25] and the Skyrme model [p6 



were, however, based on expressions for the R.H.S. of Eq. (|3.25|) which do not follow from a 
consistent identification of the isovector distribution function within the effective model. 



4 Numerical results and discussion 

The method we choose for the numerical computation of the isovector unpolarized distribu- 
tion function parallels that for computing the isosinglet unpolarized one HlOU - We evaluate 
the sums over quark single-particle levels in the background pion field by placing the soli- 
ton in a spherical 3-dimensional box, where the eigenfunctions of the Dirac Hamiltonian in 
the background pion field can be obtained by numerical diagonalization p7| . To facilitate 
the calculations we first convert the expressions Eqs.( p.l8| , |3.19|) to a spherically symmetric 
form by averaging over the 3-dimensional orientations of the spatial separation implied in 
the non-local operator Eq. (|3.1|) , see Ref.|T(| for details. In the spherically symmetric ver- 



sion of Eqs.( |3.l8| , |3.19| ) the matrix elements of the operators between single-particle quark 
states can easily be computed using the standard angular-momentum selection rules. The 
discontinuous functions of the single-particle energies and the single-particle momentum 
operator, which arise as a consequence of the presence in Eqs. (|3.18| , |3.19| ) of the functions 



5(E n + P 3 — xMn), are smoothed using the "smearing" method of Ref . [|lOj . 

For the ultraviolet regularization of distribution functions we employ the Pauli-Villars 
subtraction, Eq. (|3.21|) . For the sake of numerical stability we first evaluate the sums over 
levels Eqs.( |3.18| , |3.19p with a smooth energy cutoff, both the sums over levels obtained with 
the usual quark mass, M, and the regulator mass, Mpy. We then perform the Pauli-Villars 
subtraction, Eq.( |3.21| ), and remove the energy cutoff by extrapolation to infinity. 

We have calculated the isovector quark distribution for two values of the constituent quark 
mass, M = 350 MeV, which is the value obtained in Ref. || for the instanton vacuum, and 
M = 420 MeV. The ultraviolet cutoff is in both cases determined by fitting the pion decay 
constant, Eq.(p]), M 2 PV /M 2 = 2.52 (M = 350 MeV) and M 2 PV /M 2 = 1.90 (M = 420 MeV). 
The soliton profile, Eq.([2.3|), and the nucleon mass have been found by self-consistent min- 



imization of the Pauli-Villars regularized static energy in Ref. [11], = 1140 MeV (M = 
350 MeV) and M N = 1040 MeV (M = 420 MeV). A number of other hadronic nucleon 

(3) 

observables such as the isovector axial coupling constant g A have been calculated with this 



ultraviolet regularization in Ref.||21||. The results for the isovector unpolarized quark and 
antiquark distributions are shown in Figs.Q-^. 

Before comparing our results with the parametrizations of the experimental data it is 
instructive to study the behavior of the different contributions to the distribution function in 



3r 



The contribution Eq.(3.1S) to the isovector distribution function is a double sum over quark single- 



particle levels, similar to the moment of inertia, Eq.(2.14). In the evaluation of these quantities in a finite 



box there arises the complication that the boundary conditions which have to be imposed on the single- 
particle wave functions can lead to spurious "vacuum" contributions. In Ref.[p8[ a method has been devised 
to circumvent this problem by using two sets of basis functions in the box subject to different boundary 
conditions. We have employed this technique in the calculations of distribution functions reported here. 
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this model. Fig.[2] shows the function u(x) — d(x) for both positive and negative x, describing 
the isovector quark distribution at positive and minus the antiquark distribution at negative 
values of x. The dashed line shows the contribution of the discrete level, as defined by 
Eq.( |3.20| ). This contribution is concentrated around values x ~ 1/3 and similar in shape 
to the bound-state level contribution to the isosinglet unpolarized and isovector polarized 
distributions |TI], [II]. The contribution of the Dirac continuum, cf. Eq.( |3.2U| ), is shown by 



the dot-dashed line. Similar to the isosinglet distribution this contribution is peaked around 
x = 0; however, in the isovector case this function does not change sign at x — 0. The total 
is given by the solid line; one observes that it is essentially a smooth function except for a 
region of x around x = 0, where it is dominated by the Dirac continuum contribution. Note 
also that the calculated distribution satisfies the isospin sum rule, as discussed in Subsection 
|3~3[ (the area under the solid line in Fig.|| is unity). 

The isovector quark and antiquark distributions calculated here refer to a low normaliza- 
tion point of the order of the cutoff of the effective chiral theory (~ 600 MeV), and can be 
compared to experimental data for structure functions only after perturbative evolution to 
larger scales. Hence the "small-x" behavior of the calculated parton distribution apparent 
from Fig.^| has significance only in the sense of an input distribution at a low scale and does 
not imply a statement about the small-x behavior of the structure functions at experimental 
scales. 

The parton distribution functions calculated in the effective theory are but logarith- 
mically divergent with the UV cutoff and thus, for typical values of x, insensitive to the 
details of the UV cutoff (assuming a physically acceptable regularization scheme meeting 
the criteria discussed above). However, it is known that in the parametrically small region 



\x 



< (M/A) 2 /N c the shape of the distribution depends on the details of the UV regular- 
ization, as has recently been discussed in connection with the calculation of off-forward 
parton distributions in this approach |T3" |. For instance, it was seen there that the discon- 



tinuity at x = of the isosinglet distribution calculated with Pauli-Villars regularization 
is reduced to a smooth transition when the regularization is implemented in the form of a 
momentum-dependent constituent quark mass, as suggested by the instanton vacuum. [For 
values of x not in the vicinity of zero results are practically unchanged as compared to the 
Pauli-Villars regularization.] One may thus surmise that also the shape of the peak in the 
isovector distribution at x = depends strongly on the ultraviolet cutoff. The calculation of 
the l/iV c -suppressed distribution functions with regularization by a momentum-dependent 
mass is outside the scope of the presently available methods. The interesting general prob- 
lem of the relation of the ultraviolet cutoff of the effective theory to the "small-x" behavior 
of the low-scale parton distributions will be left for further investigation. 

In Fig.| we compare the distribution functions obtained for two different values of the 
constituent quark mass, M = 350 MeV and M = 420 MeV. With the ultraviolet regular- 
ization chosen according to Eqs. (|2.8| , p.9|) and Eq. ([3.21|) , and the soliton profile determined 



by minimization of the energy Eq. ( |2.8[ ), the only remaining free parameter in this model 
calculation is the value of the constituent quark mass, M. Rather than trying to determine 
an "optimum" value by performing a best fit to a number of hadronic observables, it is 
more interesting here to study the dependence of the calculated parton distributions on this 
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M 


level 


dx [u(x) — d 
continuum 


»] 
total 


Ig 


350 MeV 


0.030 


-0.201 


-0.171 


0.219 


420 MeV 


0.049 


-0.281 


-0.232 


0.178 



Table 1: The integral of the antiquark distribution and the Gottfried sum, Eq. ( |3.25| ), for 
constituent quark masses 350 MeV and 420 MeV. Columns 2-4: integral of the antiquark 
distribution, Jq 1 dx [u(x) —d(x)]: contribution of the discrete level, the Dirac continuum and 
the total result, cf. FigfJ. Column 5: Gottfried sum, Eq. ( |3.25| ). The experimental value 
quoted in Ref.fH is I G = 0.235 ± 0.026 at Q 2 = 4 GeV 2 . 



parameter. Fig.^| gives an idea of the variation of the results with M. 

The parton distributions computed here should be used as input for perturbative evo- 
lution, starting with a scale of the order of the cutoff (~ 600MeV). We shall not perform 
the evolution here, but rather compare our results with the parameterizations by Gliick, 
Reya, and Vogt |f[8|| . These authors generate parton distributions at experimental Q 2 from 
"valence-like" (non-singular) input quark-, antiquark and gluon distributions at a normal- 
ization point of the order of 600 MeV and obtain excellent fits to the data from deep-inelastic 
scattering and a variety of other processes. In Fig.|] we compare the isovector distribution 
of quarks, x[u(x) — d(x)}, and antiquarks, x[u(x) — d(x)], to the GRV parametrization 



One notes that the calculated distributions are systematically harder (i.e., centered at larger 
values of x) than the GRV one, indicating that the normalization point of the calculated 
distributions is even lower than that of the GRV parametrization. This is consistent with the 



conclusions drawn from the comparison of the isosinglet unpolarized || |TU , |TTJ and isovector 
polarized || [Hj distributions to the GRV parametrizations; however, in the isovector un- 
polarized case considered here the comparison of the calculated distribution with the GRV 
fit is more direct since, contrary to the isosinglet unpolarized distribution, this distribution 
does not mix with the gluon distribution under evolution, and its normalization is scale- 
independent. [Also, contrary to the isovector polarized distribution whose normalization is 

(3) 

given by the isovector axial coupling constant, g A , whose value in our approach is model- 
dependent, the normalization of the the isovector unpolarized distribution is universal thanks 
to the isospin sum rule.] 

Note that the order of magnitude and shape of the calculated isovector antiquark distri- 
bution are in good agreement with the GRV parametrization, see Fig.[|. As can be seen from 
FigPJ the antiquark distribution is dominated by the Dirac continuum contribution. We note 
that our distribution differs significantly from the result of Ref. W§1, where the contribution 



from the term Eq. ( |3.19| ) has been dropped; our distribution vanishes rapidly for large x, in 
contrast to the one of Ref.|l^], see Fig.|l]. 

With regard to the isovector antiquark distribution it is interesting to compute the Got- 
tfried sum, Eq. (|3.25|) . The results for the two constituent quark masses M = 350 MeV and 
M = 420 MeV are given in Table [I], where we also list separately the contributions of the 
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discrete level and the Dirac continuum to the integral from of the antiquark distribution. As 
can be seen, the deviation of the Gottfried sum from from 1/3 is dominated by the Dirac 



continuum contribution. The values are in reasonable agreement with the NMC value p4 



We note that, since our Dirac continuum contribution strongly differs from that of Ref. |I9 
the good agreement of the value of Iq reported there with experiment seems fortuitous. 

Finally, we note that the parton distributions studied here refer to the large- limit, 
where the nucleon is heavy, M N oc N c . The calculated distributions therefore do not go to 
zero at x = 1, rather, they are exponentially small at large x, as discussed in Ref. [[J. For this 
reason we have evaluated the integral of the antiquark distribution in Eq. ( [3.25 ) and Table [l] 
with the upper limit replaced by infinity. Numerically, the contribution from values of x > 1 
to the integral is negligible. 



5 Conclusions 

In this paper we have extended the approach to calculate parton distributions at a low 
normalization point in the large-iV c limit to the case of the distributions which appear in 
the subleading order of the l/iV c -expansion. Specifically, we have computed the isovector 
unpolarized distribution of quarks and antiquarks. The methods developed here can readily 
be generalized to compute also the other l/iV c -subleading distributions, namely the isosinglet 
longitudinally and transverse polarized ones, as well as the rotational l/iV c -corrections to 
the isovector polarized distribution. These distributions are currently being computed. 

We have found reasonable agreement of the calculated distributions with the parametriza- 
tions of the data at a low normalization point. In particular, the large- N c approach naturally 
predicts a flavor-asymmetric (i.e., non-radiative) sea quark distribution of correct sign and 
magnitude. The x-dependence of the isovector antiquark distribution compares well with 
the GRV parametrizations. It is interesting that the calculated distributions exhibit a strong 
rise at x = 0. Also, the integral of the calculated isovector antiquark distribution (the vio- 
lation of the Gottfried sum rule) is compatible with the experimental data, in view of the 
general theoretical uncertainties of the present approach. 
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A Evaluation of the isovector distribution 



In this appendix we evaluate the two contributions to the isovector distribution arising from 
the expansion in angular frequency, Eqs. fl3.llp and (|3.12| ). Our aim is to express these 
contributions in the form of sums over quark single-particle levels, which e.g. can serve as a 
starting point for a numerical calculation of the distribution function. 

We begin by rewriting the contributions Eqs. fl3.1I| ) and ( |3.12j ) in the form 

[u(x) = 

OO 

_aw e J dz0eixMNz0 J ^ j dR [^ {R) y D » {R) s*<t%*(R) 
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|0,-X) 



fi7s%R\ (A.2) 



where Tr . . . denotes the trace over Dirac and flavor indices. Here we have introduced the 
Wigner D-function in the vector representation 



D ab {R) = -TrfaRnR 1 



(A.3) 



We first have to compute the rotational matrix element. Strictly speaking, this matrix 
element contains noncommuting operators, D 3b (R) and S a , and one should be careful about 
their ordering. However, due to the average over the nucleon spin the result does not depend 
on the order, and one has 



£ / dR \&:£(R)Y D 3b (R)S a <t>T=£(R) 



E / dR [&:l(R)\ S a D 3b (R)^= s s 3 (R) 
s 3 J 



~5 ab (2T 3 ). 



(A.4) 



Using this, and passing from the time to the frequency representation for the quark Green 
functions we arrive at 



[u(x) — d(x)] 



(i) 



(2T 3 



iN c M N 



24nl 



dz v e 



0„ixM N z° 



dw 
— < 
2tt 



x Tr 



r°(l + 7 V)(z-X 



T 



lu - H{U C ) u - H(U t 



-x) 



Z 3 = -Z° ,2 I =0' 



(A.5) 
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[u(x) — d(x)] 



(2) 



-(2T 



.3, N C M N 
8tvI 



dz°e ixMNZ ° z° 



d^ c -iujz° 
2vr 



x Tr 



l+7° 7 3 )(z-X| 



u - H{U C 



-x) 



-z°, z ± =0' 



(A.6) 



[The treatment of the poles in the cu-integral implied here will be described below.] These 
expressions have the form of an integral over all space of a matrix element of the quark prop- 
agator between position eigenstates localized at different points. It is convenient to rewrite 
them as functional traces (i.e., sums over diagonal matrix elements) by introducing the fi- 
nite translation operator, which is given as the exponential of the single-particle momentum 
operator, P k , 



(z-X| 



(-X| exp(iP k z 



(A.7) 



We can then perform the integral over z° and write the result in the form 
[u(x)-d(x)] w = 

on - 

1 



, iN c Mn 



duj Sp 



r a (l + 7 V)% + ^ -xM, 



1 



AT, 



io - H(U C ) u - H(U C 



(A.i 



\u[x) 



d(x)} 



(2) 



(2T 3 ) — - — 
1 J 8nldx 



dui Sp 



1 + 7°7 3 )5(^ + P 3 - xM J 



N 



10 



H(U C ) 



(A.9) 



where Sp . . . denotes the functional trace in the space of single-particle quark states. These 
functional traces can now be computed using a basis of eigenstates of the Dirac Hamiltonian 
in the background pion field, Eq.( [2.4|) . As to the second contribution, Eq.( |A.9|) , one may 
easily show that, up to a factor, it is simply the derivative in x of the isosinglet distribution 
function, which appears in the leading order of the l/A^-expansion ||, [K|. It can be written 
simple sum over occupied single-particle levels: 



\u[x) 



d(x)} 



(2) 



-(2T 2 



,^x^[„(x) + d(x)]— 



_ (2T 3^c 



d 



{J X ^ E (n\(l + l°l 3 )6(E n + P 3 -xM N )\n). 

occup. 



(A.10) 



It is also possible to express this quantity as a sum over non-occupied levels, Eq. (|3.16|) , see 
Refs.fOl ITDfl. 
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Let us now consider the first contribution, Eq. flA.8p . It gives rise to a double sum over 
levels, 



u(x)-d(x)]V> = m) l -^f I duY,(n\T a (l+lV)^ + P 3 -xM N )\m) 

— OO 



where the sums over n, m run over all levels. We can further simplify the form of Eq. ( [A. 11 



by taking into account that the expression on the R.H.S. is real. First, this is natural, since 
the quark distribution function, u(x) — d(x), must be real. More formally, one may prove 
the reality of the R.H.S. of Eq. flA.lip starting from the locality of the Dirac Hamiltonian 
in the background field, and the fact that neither the angular velocity corrections to the 
Hamiltonian nor our ultraviolet regularization violate this property. We may thus replace 
the expression on the R.H.S. of Eq. ([A.ll|) by its real part: 

[u(x) -d(x)] {1) = 

OO 

-(2T 3 )^A Im j du V(n|r a (l + 7 7 3 )5(^ + ^ 3 -^M iV )|m) 
247r/ J 



x(m\r a \n) - (A 12) 

V ' ' ' [to - E m + iOv(E m )} [u - E n + iOv(E n )} ' V ' ; 

Here the poles in the cu-integral corresponding to the single-particle energies are shifted 
according to the occupation of the levels 

r)(E n ) = sign(£ n -£ lcv -0). (A.13) 

At E m E n we can write 
Im 



io - E m + iO<n(E m )\ [w - E n + iOr](E n )} 
-nr](E m )S(uj - E m ) 7rrj(E m )6(u - E n 



7-1 ■••|\ —III/ j-, 

io — hi n LO — Ei n 

= -7r[ V (E m )5(u - E m ) - rj(E n )5(u - E n )\ ]_ . (A.14) 

What to do in the case E m = E n depends on whether we consider the case of infinite volume, 
where most of the spectrum is continuous, or whether we work in a large but finite box, in 
which case the spectrum is quasi-discrete. We concentrate on the latter case, since it is 
relevant for the numerical calculations. In the quasi-discrete case, if E m = E n one should 
write instead of Eq. ( [A. 14 ) : 



Im 



[u - E m + iOr](E m )} [u - E n + iOrj(E n 
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1 d 
Im : — 2 = — T^Im- 



[u - E m + iQri(E m )] 2 du) u - E m + iOr](E m ) 

= 2-Tn{E m )8{w - E m ) = nv(E n )5\u-E n ). (A.15) 
ouj 

Note that the same result would be obtained if one took the limit E n — > E m in Eq.( |A.14] ). 
We can now split the double sum over m, n in Eq.( A.12 ) in two parts containing, respectively, 
terms with E m ^ E n and E m = E n and perform the integral over u using Eqs.( |A.il| ) and 
(|A.15 ). One obtains 



[u(x) — d(x)] 



(i) 



= ( 2T E F _ F (n\r a \m) 

x (m\r a (l + 7 V) {v(E n )6(E n + P 3 - xM N ) - V (E m )6(E m + P 3 - xM N )} |n) 
+ (2T 3 )^^ (n\T a \m)(m\T a (l + 1 ° 1 3 ) V (E n )5\E n + P 3 -xM N )\n). 

Ern = En 

(A.16) 

The double sum over levels in the first term here can be further simplified making use of the 
identity 

(n\T a \m)(m\r a (l + 7 ° 7 3 )<KP„ + P 3 - xM N )\n) 
= (m|r a |n)(n|r a (l + 7 7 3 )^(^n + P 3 -a;M J v)|m), (A.17) 

which holds for any pair of levels m, n. 

Equation ( A.16|) pertains to a quasi-discrete spectrum; in the continuum case one should 
replace there the summation over levels by an integral over the continuous energy, omit the 
sum over terms with E m = E n in the last line, and understand the poles (E n — Em)' 1 in the 
principal value sense. 

In Eq. ( |A.16|) the sum over n runs over all quark single-particle levels, both occupied and 
non-occupied. In order to convert it to a more standard form, particularly for use in the 
numerical calculations, we would like to rewrite the expression for the distribution function 
as a sum over either occupied or non-occupied states (cf. the standard expression for the 
moment of inertia, Eq. ( |2.14|) ). To achieve this we note that the following sum over all levels 
is zero: 

2 E F \ (n\T a \m)(m\T a (l + 7°>y 3 )6(E n + P 3 -xM N )\n) 

n,m t-'n 

Emy^En 

+ ]T (n\T a \m)(m\T a {l + 1 ° 1 3 )5\E n + P 3 -xM N )\n) = 0. (A.18) 

E rn = E n 
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To prove this identity we note that the L.H.S. can equivalently be written as the variation 
of a functional trace with a particular modified Hamiltonian 



l m J p. f dz O e i{u-xM N )J 

— oo 

x— Sp 



dA a 



r°(l + 7 V) exp(iPV) 



io - H(U C ) - AM + iO 



A=0 



dz v e 



0„-ixM N z 



o 9 



<9A a 



Sp r a (l + 7° 7 3 ) exp(iP 3 2°) exp[iz°(H(U c ) + A b r b )} 



A=0 



Due to the locality of the modified Hamiltonian, H(U C ) + A b T b , we have 

(X\exp(iP 3 z°)exp[iz°(H(U c ) +AV)]|X) = if \z 3 \ > z°. 



(A.19) 



(A.20) 



Extrapolating Eq.( |A.20| ) to z 3 — *> z° (cf. the corresponding discussion in Refs.|9|, [10[| ) we find 
that Eq.( [A.19|) is zero, which proves Eq. ( |A.18|) . Using Eq. ([A.18| ), and keeping in mind the 
definition oirj(E n ), Eq. ([A.13|) , we may now easily obtain from Eq. (|A.16|) a representation of 
the distribution function as a sum in which n runs either only over occupied states, Eq. (|3.13| ), 
or non-occupied states, Eq. (|3.15|) . We note that the equivalence of the representations of 
the distribution function as sums over occupied and non-occupied states is confirmed also 
by the numerical calculations. 

When invoking the locality condition, Eq.( |A.2(l| ), we have assumed that the ultraviolet 
regularization of the theory does not violate this property. For the regularization by Pauli- 
Villars subtraction, Eq.( p.21 ), this is indeed the case. Regularization by a cutoff, on the other 
hand, would violate Eq.( |A.20| ), which amounts to violating the anticommutation relation of 
the quark fields at space-like separations, see the discussion in Refs. 0, fTOl . 
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Figure 1: The contribution of the Dirac continuum to the isovector distribution function 
u(x) — d(x), for M = 350 MeV. At negative x the function shown describes minus the 
antiquark distribution. Solid line: Total result, as obtained from the expansion in angular 
velocity, given by the sum of Eq.( |3.18| ) and Eq. fl3.19| ). Dashed line: Result obtained dropping 
the contribution Eq. ([3.19|) , corresponding to the prescription of Ref. |19|. 
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Figure 2: Contributions to the isovector unpolarized distribution function, u(x) — d(x), for 
M = 350 MeV. At negative x the function shown describes minus the antiquark distribution. 
Dashed line: Contribution of the discrete level. Dot-dashed line: Contribution of the negative 
Dirac continuum. Solid line: Total result. 
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Figure 3: The isovector unpolarized distribution function, u(x) —d(x), for constituent quark 
masses M = 350 MeV (solid line) and M = 420 MeV (dashed line). At negative x the 
function shown describes minus the antiquark distribution. 
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Figure 4: Solid lines: The calculated isovector unpolarized quark- and antiquark distri- 
butions, x[u(x) — d(x)] and x[u(x) — d(x)], for M = 350 MeV. Shown is the total result 
corresponding to the solid line in Fig.^j. Dotted lines: The GRV NLO parametrizations 
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